EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




(cernW 



CERN-PH-EP /2010-004 

March 15, 2010 



Charged-particle multiplicities in pp interactions 
at ^fs = 900 GeV measured with the ATLAS detector 

at the LHC 



The ATLAS Collaboration 



Abstract 

The first measurements from proton-proton collisions recorded with the AT- 
LAS detector at the LHC are presented. Data were collected in December 
2009 using a minimum-bias trigger during collisions at a centre-of-mass en- 
ergy of 900 GeV. The charged-particle multiplicity, its dependence on trans- 
verse momentum and pseudorapidity, and the relationship between mean 
transverse momentum and charged-particle multiplicity are measured for 
events with at least one charged particle in the kinematic range \t]\ < 2.5 
and > 500 MeV. The measurements are compared to Monte Carlo models 
of proton-proton collisions and to results from other experiments at the same 
centre-of-mass energy. The charged-particle multiplicity per event and unit of 
pseudorapidity at 77 = is measured to be 1.333 ± 0.003 (stat.) ± 0.040 (syst.), 
which is 5-15% higher than the Monte Carlo models predict. 
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1. Introduction 



Inclusive charged-particle distributions have been measured in pp and 
pp collisions at a range of different centre-of-mass energies [l-13|. Many of 



these measurements have been used to constrain phenomenological models of 
soft-hadronic interactions and to predict properties at higher centre-of-mass 
energies. Most of the previous charged-particle multiplicity measurements 
were obtained by selecting data with a double-arm coincidence trigger, thus 
removing large fractions of diffractive events. The data were then further cor- 
rected to remove the remaining single-diffractive component. This selection 
is referred to as non-single-diffractive (NSD). In some cases, designated as 
inelastic non-diffractive, the residual double-diffractive component was also 
subtracted. The selection of NSD or inelastic non-diffractive charged-particle 
spectra involves model-dependent corrections for the diffractive components 
and for effects of the trigger selection on events with no charged particles 
within the acceptance of the detector. The measurement presented in this 
paper implements a different strategy, which uses a single-arm trigger over- 
lapping with the acceptance of the tracking volume. Results are presented 
as inclusive-inelastic distributions, with minimal model-dependence, by re- 
quiring one charged particle within the acceptance of the measurement. 

This paper reports on a measurement of primary charged particles with 
a momentum component transverse to the beam directioro > 500 MeV 
and in the pseudorapidity range \r)\ < 2.5. Primary charged particles are 
defined as charged particles with a mean lifetime r > 0.3 x 10~ 10 s directly 
produced in pp interactions or from subsequent decays of particles with a 
shorter lifetime. The distributions of tracks reconstructed in the ATLAS 
inner detector were corrected to obtain the particle-level distributions: 

1 diV ch 1 1 d 2 iV ch 1 diV ev 

and (pp) vs. n ch , 



N ev drj ' N ev 2npT dr/dpx -^" ev dn, 



ch 



where N cv is the number of events with at least one charged particle inside 
the selected kinematic range, iV c h is the total number of charged particles, 



1 The ATLAS reference system is a Cartesian right-handed co-ordinate system, with 
the nominal collision point at the origin. The anti-clockwise beam direction defines the 
positive z-axis, while the positive x-axis is defined as pointing from the collision point to 
the centre of the LHC ring and the positive y-axis points upwards. The azimuthal angle 4> 
is measured around the beam axis, and the polar angle 9 is measured with respect to the 
z-axis. The pseudorapidity is defined as r\ = — lntan(0/2). 
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n c h is the number of charged particles in an event and (pp) is the average 
Pt for a given number of charged particles. Comparisons are made to pre- 
vious measurements of charged-particle multiplicities in pp and pp collisions 
at yfs = 900 GeV centre-of-mass energies U B and to Monte Carlo (MC) 
models. 



2. The ATLAS detector 



The ATLAS detector [LjJ at the Large Hadron Collider (LHC [15J]) covers 
almost the whole solid angle around the collision point with layers of track- 
ing detectors, calorimeters and muon chambers. It has been designed to 
study a wide range of physics topics at LHC energies. For the measurements 
presented in this paper, the tracking devices and the trigger system were of 
particular importance. 

The ATLAS inner detector has full coverage in <p and covers the pseu- 
dorapidity range \rj\ < 2.5. It consists of a silicon pixel detector (Pixel), a 
silicon microstrip detector (SCT) and a transition radiation tracker (TRT). 
These detectors cover a sensitive radial distance from the interaction point 
of 50.5-150 mm, 299-560 mm and 563-1066 mm, respectively, and are im- 
mersed in a 2 T axial magnetic field. The inner-detector barrel (end-cap) 
parts consist of 3 (2x3) Pixel layers, 4 (2x9) double- layers of single-sided 
silicon microstrips with a 40 mrad stereo angle, and 73 (2x160) layers of 
TRT straws. These detectors have position resolutions of typically 10, 17 
and 130 jum for the R-<p co-ordinate and, in case of the Pixel and SCT, 115 
and 580 /im for the second measured co-ordinate. A track from a particle 
traversing the barrel detector would typically have 11 silicon hits (3 pixel 
clusters and 8 strip clusters), and more than 30 straw hits. 

The ATLAS detector has a three-level trigger system: Level 1 (LI), Level 2 
(L2) and Event Filter (EF). For this measurement, the trigger relies on the 
LI signals from the Beam Pickup Timing devices (BPTX) and the Mini- 
mum Bias Trigger Scintillators (MBTS). The BPTX are composed of beam 
pick-ups attached to the beam pipe ±175 m from the centre of the AT- 
LAS detector. The MBTS are mounted at each end of the detector in front 
of the liquid-argon end-cap calorimeter cryostats at z — ±3.56 m and are 
segmented into eight sectors in azimuth and two rings in pseudorapidity 
(2.09 < \r]\ < 2.82 and 2.82 < \rj\ < 3.84). Data were collected for this anal- 
ysis using the MBTS trigger, formed from BPTX and MBTS trigger signals. 
The MBTS trigger was configured to require one hit above threshold from 
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either side of the detector. The efficiency of this trigger was studied with a 
separate prescaled LI BPTX trigger, filtered to obtain inelastic interactions 
by inner detector requirements at L2 and EF. 



3. Monte Carlo simulation 



Low-px scattering processes may be described by lowest-order perturba- 
tive Quantum Chromodynamics (QCD) two-to-two parton scatters, where 
the divergence of the cross section at pr = is regulated by phenomeno- 
logical models. These models include multiple-parton scattering, partonic- 
matter distributions, scattering between the unresolved protons and colour 



reconnection 
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The PYTHIA 17j MC event generator implements several 
of these models. The parameters of these models have been tuned to describe 
charged-hadron production and the underlying event in pp and pp data at 
centre-of-mass energies between 200 GeV and 1.96 TeV. 

Samples of ten million MC events were produced for single-diffractive, 
double-diffractive and non- diffract ive processes using the PYTHIA 6.4.21 
generator. A specific set of optimised parameters, the ATLAS MC09 PYTHIA 
tune 18], which employs the MRST LO* parton density functions 19| and 



the pp-ordered parton shower, is the reference tune throughout this paper. 
These parameters were derived by tuning to underlying event and minimum- 
bias data from Tevatron at 630 GeVand 1.8 TeV. The MC samples generated 
with this tune were used to determine detector acceptances and efficiencies 
and to correct the data. 

For the purpose of comparing the present measurement to different phe- 
nomenological models describing minimum-bias events, the following addi- 
tional MC samples were generated: the ATLAS MC09c (lgj PYTHIA tune, 



which is an extension of the ATLAS MC09 tune optimising the strength of 
the colour reconnection to describe the (pp) distributions as a function of n c h, 
as measured by CDF in pp collisions jjij]; the PerugiaO 20] PYTHIA tune, in 
which the soft-QCD part is tuned using only minimum-bias data from the 



Tevatron and CERN pp colliders; the DW [2l| PYTHIA tune, which uses 



the virtuality-ordered showers and was derived to describe the CDF Run 
II underlying event and Drell-Yan data. Finally, the PHOJET generator 



22] was used as an alternative model. It describes low-pr physics using the 



two-component Dual Parton Model [23|, |24j, which includes soft hadronic 



processes described by Pomeron exchange and semi-hard processes described 
by perturbative parton scattering. PHOJET relies on PYTHIA for the frag- 
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mentation of partons. The versions 2 , used for this study were shown to agree 
with previous measurements (1, H, Is, 0|. 

The non-diffractive, single-diffractive and double-diffractive contributions 
in the generated samples were mixed according to the generator cross sec- 
tions to fully describe the inelastic scattering. All the events were processed 



through the ATLAS detector simulation program 25J, which is based on 



Geant4 26]. They were then reconstructed and analysed by the same pro- 
gram chain used for the data. Particular attention was devoted to the de- 
scription in the simulation of the size and position of the collision beam spot 
and of the detailed detector conditions during data taking. 

4. Event selection 

All data recorded during the stable LHC running periods between De- 
cember 6 and 15, 2009, in which the inner detector was fully operational and 
the solenoid magnet was on, were used for this analysis. During this period 
the beams were colliding head-on in ATLAS. A total of 455,593 events were 
collected from colliding proton bunches in which the MBTS trigger recorded 
one or more counters above threshold on either side. In order to perform an 
inclusive- inelastic measurement, no further requirements beyond the MBTS 
trigger and inner detector information were applied in this event selection. 
The integrated luminosity for the final event sample, which is given here for 
reference only, was estimated using a sample of events with energy deposits in 
both sides of the forward and end-cap calorimeters. The MC-based efficiency 
and the PYTHIA default cross section of 52.5 mb were then used to deter- 
mine the luminosity of the data sample to be approximately 9 /ib -1 , while 
the maximum instantaneous luminosity was approximately 5 x 10 26 cm -2 s^ 1 . 
The probability of additional interactions in the same bunch crossing was es- 
timated to be less than 0.1%. 

During this data-taking period, more than 96% of the Pixel detector, 99% 
of the SCT and 98% of the TRT were operational. Tracks were reconstructed 



offline within the full acceptance range \rj\ < 2.5 of the inner detector [271.128 



Track candidates were reconstructed by requiring seven or more silicon hits in 
total in the Pixel and SCT, and then extrapolated to include measurements 
in the TRT. Typically, 88% of tracks inside the TRT acceptance (\r)\ < 



2 PHOJET 1.12 with PYTHIA 6.4.21 
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2) include a TRT extension, which significantly improves the momentum 
resolution. 

This paper reports results for charged particles with p^ > 500 MeV, which 
are less prone than lower-px particles to large inefficiencies and their associ- 
ated systematic uncertainties resulting from interactions with material inside 
the tracking volume. To reduce the contribution from background events and 
non-primary tracks, as well as to minimise the systematic uncertainties, the 
following criteria were required: 



• the presence of a primary vertex [29[ reconstructed using at least three 
tracks, each with: 

— p T > 150 MeV, 

— a transverse distance of closest approach with respect to the beam- 
spot position \d>Q S \ < 4 mm. 

• at least one track with: 

— p T > 500 MeV, 

— a minimum of one Pixel and six SCT hits, 

— transverse and longitudinal impact parameters calculated with 
respect to the event primary vertex \do\ < 1.5 mm and \zq\ ■ 
sin 8 < 1.5 mm, respectively. 

These latter tracks were used to produce the corrected distributions and 
will be referred to as selected tracks. The multiplicity of selected tracks 
within an event is denoted by ns e i- In total 326,201 events were kept after 
this offline selection, which contained 1,863,622 selected tracks. The inner 
detector performance is illustrated in Fig. [T] using selected tracks and their 
MC simulation. The shapes from overlapping Pixel and SCT modules in the 
forward region and the inefficiency from a small number of disabled Pixel 
modules in the central region are well modelled by the simulation. The 
simulated impact-parameter distributions describe the data to better than 
10%, including their tails as shown in the inserts of Fig. [Tb and d. The 
difference between data and MC observed in the central region of the do 
distribution is due to small residual misalignments not simulated in the MC, 
which are found to be unimportant for this analysis. 
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Figure 1: Comparison between data (dots) and minimum-bias ATLAS MC09 simulation 
(histograms) for the average number of Pixel hits (a) and SCT hits (b) per track as a 
function of ry, and the transverse (c) and longitudinal (d) impact parameter distributions 
of the reconstructed tracks. The MC distributions in (c) and (d) are normalised to the 
number of tracks in the data. The inserts in the lower panels show the distributions in 
logarithmic scale. 

Trigger and vertex-reconstruction efficiencies were parameterized as a 
function of the number of tracks passing all of the track selection require- 
ments except for the constraints with respect to the primary vertex. Instead, 
the transverse impact parameter with respect to the beam spot was required 
to be less than 4 mm, which is the same requirement as that used in the 
primary vertex reconstruction preselection. The multiplicity of these tracks 
in an event is denoted by nfjj. 

5. Background contribution 

There are two possible sources of background events that can contaminate 
the selected sample: cosmic rays and beam-induced background. A limit on 
the fraction of cosmic-ray events recorded by the LI MBTS trigger during 
data taking was determined from cosmic-ray studies, the maximum number 



7 



of proton bunches, and the central trigger processor clock width of 25 ns, and 
was found to be smaller than 10~ 6 . Beam-induced background events can 
be produced by proton collisions with upstream collimators or with residual 
particles inside the beam pipe. The LI MBTS trigger was used to select 
beam-induced background events from un-paired proton bunch-crossings. By 
applying the analysis selection criteria to these events, an upper limit of 10 -4 
was determined for the fraction of beam-induced background events within 
the selected sample. The requirement of a reconstructed primary vertex is 
particularly useful to suppress the beam-induced background. 

Primary charged-particle multiplicities are measured from selected-track 
distributions after correcting for the fraction of secondary particles in the 
sample. The potential background from fake tracks is found to be less 
than 0.1% from simulation studies. Non-primary tracks are mostly due 
to hadronic interactions, photon conversions and decays of long-lived par- 
ticles. Their contribution was estimated using the MC prediction of the 
shape of the do distribution, after normalising to data for tracks within 
2 mm < | do | < 10 mm, i.e. outside the range used for selecting tracks. 
The simulation was found to reproduce the tails of the d distribution of the 
data, as shown in Fig. [Tb, and the normalisation factor between data and 
MC was measured to be 1.00 ±0.02 (stat.) ±0.05 (syst.). The MC was then 
used to estimate the fraction of secondaries in the selected-track sample to 
be (2.20 ±0.05 (stat.) ±0.11 (syst.))%. This fraction is independent of ns e i, 
but shows a dependence on pi* and a small dependence on rj. While the 
correction for secondaries was applied in bins of p^, the dependence on r\ was 
incorporated into the systematic uncertainty. 

6. Selection efficiency 

The data were corrected to obtain inclusive spectra for charged pri- 
mary particles satisfying the event-level requirement of at least one primary 
charged particle within px > 500 MeV and \rj\ < 2.5. These corrections in- 
clude inefficiencies due to trigger selection, vertex and track reconstruction. 
They also account for effects due to the momentum scale and resolution, and 
for the residual background from secondary tracks. 

Trigger efficiency. The trigger efficiency was measured from an independent 
data sample selected using the control trigger introduced in Section 2. This 
control trigger required more than 6 Pixel clusters and 6 SCT hits at L2, 
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and one or more reconstructed tracks with > 200 MeV at the EF. The 
vertex requirement for selected tracks was removed for this study, to avoid 
correlations between the trigger and vertex- reconstruction efficiencies for LI 
MBTS triggered events. The trigger efficiency was determined by taking the 
ratio of events from the control trigger in which the LI MBTS also accepted 
the event, over the total number of events in the control sample. The result 
is shown in Fig. [2^ as a function of ngjj. The trigger efficiency is nearly 100% 
everywhere and the requirement of this trigger does not affect the px and r\ 
track distributions of the selected events. 
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Figure 2: Trigger (a) and vertex-reconstruction (b) efficiencies as a function of the variable 
riggj defined in Section^ track-reconstruction efficiency as a function of px (c) and of 77 (d). 
The vertical bars represent the statistical uncertainty, while the shaded areas represent 
the statistical and systematic uncertainties added in quadrature. The two bottom panels 
were derived from the PYTHIA ATLAS MC09 sample. 



Vertex-reconstruction efficiency. The vertex-reconstruction efficiency was de- 
termined from the data, by taking the ratio of triggered events with a re- 
constructed vertex to the total number of triggered events. It is shown in 
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Fig. U3 as a function of nfjj. The efficiency amounts to approximately 67% 
for the lowest bin and rapidly rises to 100% with higher multiplicities. The 
dependence of the vertex-reconstruction efficiency on the n and p T of the se- 
lected tracks was studied. The n dependence was found to be approximately 
flat for Tigjj > 1 and to decrease at larger n for events with n^\ = 1. This 
dependence was corrected for. No dependence on p^ was observed. 

Track-reconstruction efficiency. The track-reconstruction efficiency in each 
bin of the Pt-tj acceptance was determined from MC. The comparison of 
the MC and data distributions shown in Fig. [T] highlights their agreement. 
The track-reconstruction efficiency was defined as: 

, , iV r r chcd (PT,V) 

CbintPT,//) = , r 7 v— , 

N g en{PT,V) 

where pt and n are generated quantities, and iV™ ) ? tched (pr ) v) an d N gen (j)T, v) 
are the number of reconstructed tracks in a given bin matched to a generated 
charged particle and the number of generated charged particles in that bin, 
respectively. The matching between a generated particle and a reconstructed 
track was done using a cone-matching algorithm in the rj-(j) plane, associ- 
ating the particle to the track with the smallest AR = a/ (A0) 2 + (An) 2 
within a cone of radius 0.05. The resulting reconstruction efficiency as a 
function of px integrated over n is shown in Fig. [2b. The drop to ~ 70% for 
Pt < 600 MeV is an artefact of the pr cut at the pattern-recognition level 
and is discussed in Section |HJ The reduced track-reconstruction efficiency in 
the region \rj\ > 1 (Fig. [2U) is mainly due to the presence of more material in 
this region. These inefficiencies include a 5% loss due to the track selection 
used in this analysis, approximately half of which is due to the silicon-hit 
requirements and half to the impact-parameter requirements. 



7. Correction procedure 

The effect of events lost due to the trigger and vertex requirements can 
be corrected for using an event-by-event weight: 

fc tngl"'Sely | fcytx^SelJ 

where etri g (^ B c S i) an d e vtx(^ B o S i) are ^ ne trigger and vertex reconstruction effi- 
ciencies discussed in Section [6j The vertex- reconstruction efficiency for events 



10 



with rigjj = 1 includes an 77- dependent correction which was derived from the 
data. 

The and rj distributions of selected tracks were corrected on a track- 
by-track basis using the weight: 

W tT k(PT,V) = r T ' ( X ~ /scc(Pt)) ■ (1 - /okr(PT,?7)), 

where e^n is the track-reconstruction efficiency described in Section O and 
f S ec{.PT) is the fraction of secondaries determined as described in Sectional 
The fraction of selected tracks for which the corresponding primary particles 
are outside the kinematic range, / kr(pr, v)i originates from resolution effects 
and has been estimated from MC. Bin migrations were found to be due solely 
to reconstructed track momentum resolution and were corrected by using the 
resolution function taken from MC. 

In the case of the distributions versus ri c h, a track-level correction was 



applied by using Bayesian unfolding 30] to correct back to the number of 



charged particles. A matrix M^seh which expresses the probability that 
a multiplicity of selected tracks ns c i is due to n cil particles, was populated 
using MC and applied to obtain the n c h distribution from the data. The 
resulting distribution was then used to re-populate the matrix and the cor- 
rection was re-applied. This procedure was repeated without a regularisa- 
tion term and converged after four iterations, when the change in the dis- 
tribution between iterations was found to be less than 1%. It should be 
noted that the matrix cannot correct for events which are lost due to track- 
reconstruction inefficiency. To correct for these missing events, a correction 
factor 1/(1 — (1 — e(n c h))" ch ) was applied, where e(n c h) is the average track- 
reconstruction efficiency. 

In the case of the (pt) versus distribution, each event was weighted 
by w cv(nici)- For each ns e i a MC-based correction was applied to convert 
the reconstructed average pr to the average of primary charged particles. 
Then the matrix M c h,sei was applied as described above. 



8. Systematic uncertainties 

Numerous detailed studies have been performed to understand possible 
sources of systematic uncertainties. The main contributions are discussed 
below. 
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Trigger. The trigger selection dependence on the pr and rj distributions of 
reconstructed tracks was found to be flat within the statistical uncertainties 
of the data recorded with the control trigger. The statistical uncertainty 
on this result was taken as a systematic uncertainty of 0.1% on the overall 
trigger efficiency 

Since there is no vertex requirement in the data sample used to measure 
the trigger efficiency, it is not possible to make the same impact-parameter 
cuts as are made on the final selected tracks. Therefore the trigger efficiency 
was measured using impact-parameter constraints with respect to the pri- 
mary vertex or the beam spot and compared to that obtained without such 
a requirement. The difference was taken as a systematic uncertainty of 0.1% 
for nfj < 3. 

The correlation of the MBTS trigger with the control trigger used to select 
the data sample for the trigger-efficiency determination was studied using the 
simulation. The resulting systematic uncertainty was found to affect only the 
case rigjj = 1 and amounts to 0.2%. 

Vertex reconstruction. The run-to-run variation of the vertex-reconstruction 
efficiency was found to be within the statistical uncertainty. The contribu- 
tion of beam-related backgrounds to the sample selected without a vertex 
requirement was estimated by using non-colliding bunches. It was found to 
be 0.3% for 

n Sci = 1 an d smaller than 0.1% for higher multiplicities, and was 
assigned as a systematic uncertainty. This background contribution is larger 
than that given in Section |5l since a reconstructed primary vertex was not 
required for these events. 

Track reconstruction and selection. Since the track-reconstruction efficiency 
is determined from MC, the main systematic uncertainty results from the 
level of disagreement between data and MC. 

Three different techniques to associate generated particles to reconstructed 
tracks were studied: a cone-matching algorithm, an evaluation of the fraction 
of simulated hits associated to a reconstructed track and an inclusive tech- 
nique using a correction for secondary particles. A systematic uncertainty of 
0.5% was assigned from the difference between the cone-matching and the 
hit-association methods. 

A detailed comparison of track properties in data and simulation was per- 
formed by varying the track-selection criteria. The largest deviations between 
data and MC were observed by varying the zq ■ sin 8 selection requirement, 
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and by varying the constraint on the number of SCT hits. These deviations 
are generally smaller than 1% and rise to 3% at the edges of the 77 range. 

The systematic effects of misalignment were studied by smearing simu- 
lation samples by the expected residual misalignment and by comparing the 
performance of two alignment algorithms on tracks reconstructed from the 
data. Under these conditions the number of reconstructed tracks was mea- 
sured and the systematic uncertainty on the track reconstruction efficiency 
due to the residual misalignment was estimated to be less than 1%. 

To test the influence of an imperfect description of the detector material 
in the simulation, two additional MC samples with approximately 10% and 
20% increase in radiation lengths of the material in the Pixel and SCT active 
volume were used. The impact of excess material in the tracking detectors 
was studied using the tails of the impact-parameter distribution, the length 
of tracks, and the change in the reconstructed K$ mass as a function of the 
decay radius, the direction and the momentum of the K$. The MC with 
nominal material was found to describe the data best. The data were found 
to be consistent with a 10% material increase in some regions, whereas the 
20% increase was excluded in all cases. The efficiency of matching full tracks 
to track segments reconstructed in the Pixel detector was also studied. The 
comparison between data and simulation was found to have good agreement 
across most of the kinematic range. Some discrepancies found for \r]\ > 1.6 
were included in the systematic uncertainties. From all these studies a sys- 
tematic uncertainty on the track reconstruction efficiency of 3.7%, 5.5% and 
8% was assigned to the pseudorapidity regions \r]\ < 1.6, 1.6 < \rj\ < 2.3 and 
\r]\ > 2.3, respectively. 

The track reconstruction efficiency shown in Fig. [2b rises sharply in the 
region 500 < p-r < 600 MeV. The observed turn-on curve is produced by 
the initial pattern recognition step of track reconstruction and its associated 
Pt resolution, which is considerably worse than the final resolution. The 
consequence is that some particles which are simulated with p^ > 500 MeV 
are reconstructed with momenta below the selection requirement. This effect 
reduces the number of selected tracks. The shape of the threshold was studied 
in data and simulation and a systematic uncertainty of 5% was assigned to 
the first px bin. 

In conclusion, an overall relative systematic uncertainty of 4.0% was as- 
signed to the track reconstruction efficiency for most of the kinematic range 
of this measurement, while 8.5% and 6.9% were assigned to the highest |?7| 
and to the lowest bins, respectively. 
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Momentum scale and resolution. To obtain corrected distributions of charged 
particles, the scale and resolution uncertainties in the reconstructed px and 
r) of the selected tracks have to be taken into account. Whereas the uncer- 
tainties for the rj measurement were found to be negligible, those for the px 
measurement are in general more important. The inner detector momentum 
resolution was taken from MC as a function of px and 77. It was found to vary 
between 1.5% and 5% in the range relevant to this analysis. The uncertainty 
was estimated by comparing with MC samples with a uniform scaling of 10% 
additional material at low px and with large misalignments at higher px- 
Studies of the width of the mass peak for reconstructed K$ candidates in the 
data show that these assumptions are conservative. The reconstructed mo- 
mentum scale was checked by comparing the measured value of the K§ mass 
to the MC. The systematic uncertainties from both the momentum resolution 
and scale were found to have no significant effect on the final results. 

Fraction of secondaries. The fraction of secondaries was determined as dis- 
cussed in Section The associated systematic uncertainty was estimated 
by varying the range of the impact parameter distribution that was used to 
normalise the MC, and by fitting separate distributions for weak decays and 
material interactions. The systematic uncertainty includes a small contri- 
bution due to the rj dependence of this correction. The total uncertainty is 
0.1%. 

Correction procedure. Several independent tests were made to assess the 
model dependence of the correction matrix M c ^e\ and the resulting sys- 
tematic uncertainty. In order to determine the sensitivity to the px and rj 
distributions, the matrix was re-populated using the other MC parameteriza- 
tions described in Section [3] and by varying the track-reconstruction efficiency 
by ± 5%. The correction factor for events lost due to the track-reconstruction 
inefficiency was varied by the same amount and treated as fully correlated. 
For the overall normalisation, this leads to an uncertainty of 0.4% due to 
the model dependence and of 1.1% due to the track-reconstruction efficiency. 
The size of the systematic uncertainties on n c h increases with the multiplicity. 

The correction for the (px) was also studied using the different PYTHIA 
tunes and PHO JET. The change was found to be less than 2% over the whole 
sample. 

As the track-reconstruction efficiency depends on the particle type, the 
uncertainty in the composition of the charged particles in the minimum-bias 
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MC sample was studied. The relative yields of pions, kaons and protons in 
the simulation were separately varied by ±10%. These variations, combined 
with changing the fraction of electrons and muons by a factor of three, re- 
sulted in a systematic uncertainty of 0.2%. 

The systematic uncertainty on the normalisation and on the number of 
charged particles were treated separately. In each of these two groups the sys- 
tematic uncertainties were added in quadrature. These were then combined 
taking into account their anti-correlation and were propagated to the final 
distributions. Table [TJ summarises the various contributions to the systematic 
uncertainties on the charged-particle density at rj = 0. 



Systematic uncertainty on the number of events, iV ev 


Trigger efficiency 




< 0.1% 


Vertex-reconstruction efficiency 




< 0.1% 


Track-reconstruction efficiency 




1.1% 


Different MC tunes 




0.4% 


Total uncertainty on N ev 


1.2% 


Systematic uncertainty on (l/A^ ev )"(dA^ c h 


/drj) at rj — 


Track-reconstruction efficiency 




4.0% 


Trigger and vertex efficiency 




< 0.1% 


Secondary fraction 




0.1% 


Total uncertainty on 7V ev 




-1.2% 


Total uncertainty on (l/N ev )-(dN ch /dr]) at r] = 





2.8% 



Table 1: Summary of systematic uncertainties on the number of events, N cv , and on 
the charged-particle density (l/iV ov )-(dA^ c h/d77) at i] = 0. The uncertainty on A^ cv is 
anticorrelated with dN c \i/dr). All other sources are assumed to be uncorrelated. 



9. Results 

The corrected distributions for primary charged particles for events with 
^ch > 1 in the kinematic range pr > 500 MeV and \rj\ < 2.5 are shown in 
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Fig. |3l where they are compared to predictions of models tuned to a wide 
range of measurements. The data are presented as inclusive-inelastic distri- 
butions with minimal mo del- dependent corrections to facilitate the compar- 
ison with models. 

The charged-particle pseudorapidity density is shown in Fig. |3^. It is 
measured to be approximately flat in the range \r)\ < 1.5, with an average 
value of 1.333 ± 0.003 (stat.) ±0.040 (syst.) charged particles per event and 
unit of pseudorapidity in the range \rj\ < 0.2. The particle density is found 
to drop at higher values of \rj\. All MC tunes discussed in this paper are 
lower than the data by 5-15%, corresponding to approximately 1-4 standard 
deviations. The shapes of the models are approximately consistent with the 
data with the exception of PYTHIA DW. 

The iV c h distribution in bins of pr is shown in Fig. [3b and is constructed 
by weighting each entry by 1/pr- The MC models do not reproduce the data 
for pt > 0.7 GeV. The most significant difference is seen for the PHOJET 
generator. 

The multiplicity distribution as a function of n c ^ is shown in Fig. [3b. The 
PYTHIA models show an excess of events with n c ^ = 1 with respect to the 
data, while the fraction of events with n c h > 10 is consistently lower than 
the data. The net effect is that the integral number of charged particles 
predicted by the models are below that of the data (Fig. [3^ and [3b). The 
PHOJET generator successfully models the number of events with n c h = 1, 
while it deviates from the data distributions at higher values of n c h- 

The average p^ as a function of n c h is illustrated in Fig. [3bl. It is found 
to increase with increasing n c h and a change of slope is observed around 
n c h = 10. This behaviour was already observed by the CDF experiment 
in pp collisions at 1.96 TeV The PerugiaO parameterization, which was 
tuned using CDF minimum-bias data at 1.96 TeV, describes the data well. 
The other models fail to describe the data below n c h ~ 25, with the exception 
of the PYTHIA-MC09c tune. 

The iV c h distribution as a function of p^ in the kinematic range p^ > 
500 MeV and \r]\ < 2.5 is shown in Fig. HI The CMS [l| results at the same 
centre-of-mass energy are superimposed. The number of charged particles 
in the CMS data is consistently lower than the data presented in this pa- 
per. This offset is expected from the CMS measurement definition of NSD 
events, where events with = enter the normalisation and the number 
of lower transverse momentum particles are reduced by the subtraction of 
the PYTHIA single diffractive component. The UAl [5] results, normalised 
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Figure 3: Charged-particle multiplicities for events with n c h > 1 within the kinematic 
range pt > 500 MeV and \rj\ < 2.5. The panels show the charged-particle multiplicity as a 
function of pseudorapidity (a) and of the transverse momentum (b), the charged-particle 
multiplicity (c), and the average transverse momentum as a function of the number of 
charged particles in the event (d). The dots represent the data and the curves the predic- 
tions from different MC models. The vertical bars represent the statistical uncertainties, 
while the shaded areas show statistical and systematic uncertainties added in quadrature. 
The values of the ratio histograms refer to the bin centroids. 
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by their associated cross section measurement, are also overlaid. They are 
approximately 20% higher than the present data. A shift in this direction is 
expected from the double-arm scintillator trigger requirement used to collect 
the UAl data, which rejected events with low charged-particle multiplicities. 

To compare more directly the present data with results from CMS, the 
mean charged-particle density was calculated in the range \rj\ < 2.4 and a 
model dependent correction was applied to form an NSD particle density. 
For the calculation of the NSD value the PYTHIA DW tune was selected 
due to its similarity with the tune used in the CMS analysis. This genera- 
tor set-up was used to produce a correction for the removal of the fraction 
of single diffractive events, the removal of electrons from tt° Dalitz decays 
and the addition of non-single diffractive events with no charged particles 
within the kinematic range pr > 500 MeV and \rj\ < 2.5. The net effect 
of the correction is to reduce the charged-particle multiplicity. The result- 
ing value 1.240 ± 0.040 (syst.) is consistent with the CMS measurement of 

I. 202 ± 0.043 (syst.) in the kinematic range of p T > 500 MeV and \r)\ < 2.4. 

10. Conclusions 

Charged-particle multiplicity measurements with the ATLAS detector 
using the first collisions delivered by the LHC during 2009 are presented. 
Based on over three hundred thousand proton-proton inelastic interactions, 
the properties of events with at least one primary charged particle produced 
within the kinematic range \r]\ < 2.5 and > 500 MeV were studied. 
The data were corrected with minimal model dependence to obtain inclu- 
sive distributions. The charged-particle multiplicity per event and unit of 
pseudorapidity at rj = is measured to be 1.333±0.003(stat.)±0.040(syst.), 
which is 5-15% higher than the Monte Carlo model predictions. The selected 
kinematic range and the precision of this analysis highlight clear differences 
between Monte Carlo models and the measured distributions. 
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